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Abstract 

We have shown elsewhere that the presence of mixed-culture growth of microbial species in fermentation processes 
can be detected with high accuracy by employing the wavelet transform. This is achieved because the crosses in 
the different growth processes contributing to the total biomass signal appear as singularities that are very well 
evidenced through their singularity cones in the wavelet transform. However, we used very simple two-species 
cases. In this work, we extend the wavelet method to a more complicated illustrative fermentation case of three 
microbial species for which we employ several wavelets of different number of vanishing moments in order to 
eliminate possible numerical artifacts. Working in this way allows to filter in a more precise way the numerical 
values of the Holder exponents. Therefore, we were able to determine the characteristic Holder exponents for the 
corresponding crossing singularities of the microbial growth processes and their stability logarithmic scale ranges 
up to the first decimal in the value of the characteristic exponents. Since calibrating the mixed microbial growth 
by means of their Holder exponents could have potential industrial applications, the dependence of the Holder 
exponents on the kinetic and physical parameters of the growth models remains as a future experimental task. 
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Resumen 

Hemos mostrado en un trabajo anterior que la presencia de crecimiento mixto de poblaciones en procesos de 
fermentacion puede ser inferida por medio del analisis ondeleta, el cual detecta con mucha precision cruces en la 
serial de biomasa total. Esto fue realizado debido a que los cruces en los diferentes procesos de crecimiento que 
contribuyen en la serial de biomasa total surjan como singularidades muy bien delineadas en sus conos de singu- 
laridad de la transformada ondeleta. Sin embargo, usamos casos muy simples de dos especies. En este trabajo 
extendemos el estudio de tipo ondeleta a un caso ilustrativo de fermentacion mas complicado de tres poblaciones 
microbianas, en el cual se usaron varios tipos de ondeletas de diferente niimero de momentos de desvanecimiento 
para evitar los posibles artefactos numericos. Trabajando de esta forma nos permite filtrar de manera mas precisa 
los valores numericos de los exponentes de Holder. Por consiguiente, pudimos determinar los exponentes carac- 
terfsticos de Holder para las singularidades de cruce correspondientes a los procesos de crecimiento microbiano 
asi como de sus ranges estables del logaritmo de la escala hasta el primer decimal en el valor de los exponentes 
caracterfsticos. Debido a que la calibracion del crecimiento mixto microbiano por medio de los exponentes de 
Holder podria tener aplicaciones industriales potenciales, la dependencia de los exponentes de Holder en funcion 
de los parametros cineticos de los modelos de crecimiento queda como una futura tarea experimental. 

Palabras claves: crecimiento de cultivo mixto, biomasa total, ondeletas, exponentes de Holder. 
Archivo fuente: jart2010.tex 

1 Introduction 

The growth of microbial species in media containing two or more substrates hmiting the growth is of consid- 
erable biotechnological and bioengineering interest [1] . The mixed growth of microorganisms occurs in many 
industrial processes. An important class of such processes is that of the traditional fermentation food and 
drinks in which either ambient microorganisms or inoculums with selected microorganisms are used. The 
presence of different microbial species and substrates is a dominant factor for the quality and quantity of 
the final product. Biomass measurements are among the most important measurements in any cell culture 
process and detecting the presence or absence of the mixed culture growing processes is a demanding issue. 
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In general terms, the biomass measurements can by divided in continuous and on-line, and discrete and 
off-line. The continuous biomass acquisition signals and reliable on-line methods yield valuable knowledge 
on the status of the process and can facilitate the process monitoring and control. Several methods have 
been developed for the continuous measurement of biomass. Traditional optical methods such as absorbance 
measurements are probably the easiest ways and are the most commonly used. Additionally, capacitance 
has been also used to estimate the biomass concentration in a continuous and online way, see for instance 
the web pages of the Fogale and Applikon companies: http://www.fogale.fr/biotech/pages/yeasts-2.php and 
http://www.applikonbio.com/applikonbio/c8-l.htm, respectively. Our main goal is to show that it is possi- 
ble to infer the mixed culture microbial growth based only on the knowledge of total biomass data, without 
employing complicated techniques. The alternative procedure that we emphasized recently and it is here 
shortly presented is to apply the wavelet method to the fermentation data with the purpose to detect sin- 
gularities in the growth curves (signals). In this case, one can consider the mixed culture growth curves as 
more or less regular signals containing crosses due to the difference in growth of the different species. In the 
wavelet literature, one can find fundamental papers in which it is shown that the wavelet techniques provide 
an efficient tool for detecting singularities, see for example Mallat and Hwang [2]. 



2 Microbial kinetics 

The fermentation processes on which we focus here evolve in stirred batch tank reactors. In such a case, the 
biomass converts the substrate in additional biomass and products. The dynamical model of the bioreactor 
for biomass production is given by the following ordinary differential equations: 

dxiA 



dt 



XlA ^J.i {X2,i) 



dX2,i Xit 

xi_i are the biomass concentrations, X2,i the concentrations of the substrate, Yi is the yield factor defined as 
the ratio between the amount of produced biomass per unit of time and the amount of consumed substrate 
per unit of time for each species, fJ-i{x2,i) are the growth rates of each species and the subindex i represents 
the number of species and substrates involved in the fermentation process. The growth rates relate the 
changes in the biomass concentrations with the (negative) changes of the substrate concentrations per unit 
of time. We will use the formulas of Monod and Haldane, respectively: 

1. - For the substrate saturation model (the Monod equation): 

/ \ ^maxi X2 i 

K-i, + X2,i 

where Ki. is the saturation (or Monod) constant. 

2. - For the substrate inhibition model (the Haldane equation): 



Kl, + X2A + K2,xli 

where in addition to the saturation constant there is an inhibition constant K2i ■ 

In both models, fJ-maxi is the maximal specific growth rate. For simplicity, this parameter has been taken 
as unity in the illustrative case presented in the following. The value of Ki. expresses the affinity of the 
biomass for the substrate. The Monod growth kinetics can be considered as a special case of the substrate 
inhibition kinetics for if 2- = 0, i.e., when the inhibition terms vanish. The total biomass signal (TBS) is 
naturally given by 

m 

2/W= (1) 

i=l 

where m is the total number of microbial species inside the bioreactor. In our previous work [7, we considered 
the simple case of two species and two substrates. In Section 3, we will present a more complicated case, 
namely the case of three species, i.e., to = 3, with identical initial conditions for all species and substrates. 
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3 Wavelet analysis 



In the context of signal theory the total biomass signal y{t) is projected onto a wavelet basis through the 
correlation of the signal with an integral kernel which is a scaled and translated version of an analyzing 
wavelet 

1 C°° ft-b\ 
T4y]{a,b)^^J y{t)i^i^-^jdt, (2) 

where ip is the analyzing wavelet. The values of Eq. 1^ are called the wavelet coefficients of the signal y(t). 
In general, a signal can be disentangled into three parts: (a) a mostly smooth continuous structure that 
can be represented by low-order (piecewise) polynomials, (b) a discontinuous part intrinsic to the particular 
nature of the signal, and (c) a noise part which can also be quite discontinuous. In general, for any function 
f{t) we can define its moments as Mk = t'' f{t)dt, for k a nonnegative integer. In the case of wavelets, 
an analyzing wavelet ijj has p vanishing moments if and only if it satisfies 

/oo 
t''ij(t)dt = , for aU < fc < p . (3) 
-oo 

This property means that the analyzing wavelet ijj is orthogonal to any polynomial having degree up to 
p — 1. The closer a signal is to being orthogonal to the analyzing wavelet, the smaller the resulting wavelet 
coefficient. Thus, if the analyzed signal can be approximated by a low-order polynomial over the support 
of the wavelet, then the wavelet is nearly orthogonal to the signal and the resulting wavelet coefficients 
will be very small (resulting from the approximation error). The smaller the wavelet support, the greater 
the precision in detecting discontinuities and the shorter the processing time. Moreover, the support of a 
wavelet is directly proportional to its number of vanishing moments. We will display results presenting the 
same total biomass signal analyzed with wavelets of the Gaussian family for different vanishing moments, 
from two to four. In addition, as in our previous work we use the wavelet transform modulus maxima 
procedure introduced by Mallat and Hwang [2] that provides a precise location of the singularity and a 
clear cut measure of its Holder exponent (HE) from the scaling of the wavelet transforms along the so-called 
modulus maxima lines on which the transforms reach local maxima with respect to the position coordinate. 
For more details, the reader is directed to the works |2j and [3]. 



4 A case study: Two Haldane and one Monod species 

Depending on the saturation and inhibition parameters and the types of fermentative growth, the crosses in 
the case of three growing species can occur at well-separated moments. For example, the total biomass signal 
may have two separated cones of influence in its wavelet transform. Such a case for two Haldane populations 
and one of Monod type is displayed in the figures of this work for different analyzing wavelets. Processing 
the same signal with several analyzing wavelets is recommendable when a better identification of the HEs 
is sought and helps to get rid of all sorts of artifacts. We used the second, third, and fourth derivatives 
of the Gaussian function as analyzing wavelets. It is well known that these Gaussian derivative wavelets 
have the same number of vanishing moments as the order of the derivative. In each case the range of scales 
with the same HE up to the first decimal digit in all three continuous wavelet transforms have been found. 
The same type of analysis, namely using wavelets of different vanishing moments to sense a given signal, 
have been used by Arneodo et o/ [3] in a research of the fractal properties of DNA genome sequences. The 
argument for working in this way is that the common features showing up in all wavelet transforms belong 
with much higher confidence to the signal itself and it is not a numerical artifact caused by the integral 
transform. Following this idea, in the case of three microbial populations presented here, we can conclude 
that the HE is 1.6 for the first crossing singularity and 1.3 for the second one. This is so because these values 
are preserved in sufficiently long scale ranges although not for the entire scale interval that we considered. 

• The first crossing singularity 

The results concerning the first encountered singularity are displayed in the figures 1,2, and 3. From 
the plots (c) in each of the figures, we notice that we are able to get a stable first decimal digit of the 
HE if the low scales in the calculation of the wavelet transform are eliminated roughly proportional 
with the order of the Gaussian derivative employed. 
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Fig. 1: (a) Total biomass signal (TBS) for three microbial species with different growth rates: two are of Haldane 
type with Ki^ = 0.02 g/1 and — 0.04 g/1, and Ki^ = 0.03 g/1 and K22 ~ 0.9 g/1, respectively and the third 
one is of Monod type with Ki = 0.09g/l. (b) The singularity analysis has been performed with the second Gaussian 
derivative wavelet, (c) The HE a is shown in two different scale ranges. 




(a) 
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Fig. 3: The same case as in previous figures but with the singularity analysis performed with the fourth Gaussian 
derivative wavelet. 



• The second crossing singularity 

The results concerning the second singularity are displayed in table 1 (we do not include the plots for 
this singularity here). In this case, we could not assess a similar trend for the low scales. However, we 
identified the scale ranges on which we got a stable first decimal digit of the exponent for the three 
Gaussian derivatives, as one can see from the right plots (c). 

• The plateau singularity 

The plateau singularity is less interesting because it indicates only the end of the fermentation process. 
The only useful fact is its detection and we noticed that it is indeed detected with any of the employed 
wavelets. 

We collect the results of our illustrative example in table 1. One can notice that we can determine a 
range of scales where the HE is the same up to the first decimal digit independently of the number of 
vanishing moments of the used wavelet. In the case of the first crossing, HE is constant only over two 
log scales (5-6), whereas for the second crossing it stays constant over three log scales (3-5). 



Type of wavelet 


1st crossing 


2nd crossing 


Plateau 


Gaussian derivatives 


Rif 


i?2/ 


-Rls R2s 


Detected ? 


G2 


1.5220 


1.6028 


1.1600 


1.3200 


Yes 


G3 


1.8100 


1.6357 


1.2440 


1.3730 


Yes 


G4 


2.2320 


1.6533 


1.2429 


1.3514 


Yes 



Table 1: The Holder exponents obtained in two different ranges of scales Ri and R2 for the first (subindex /) and 
the second (subindex s) crossing singularities, respectively, with the wavelet functions corresponding to the second 
(G2), third (G3) and fourth (G4) derivatives of the Gaussian function. 



As a general remark, notice that we did not use the first Gaussian derivative as an analyzing wavelet 
because in this case the HEs are very close to unity which means that the signal is continuous at the 
singularity detected with higher-order derivative Gaussian wavelets; this fact tells us that the continu- 
ous wavelet transform with this wavelet cannot detect the singularities in the signal independently of 
the range of scales that one takes into account, see also Mallat and Hwang [2j. 
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5 Conclusions 



The presence of mixed culture growth in fermentation processes can be inferred starting from the data of 
total biomass processed with the wavelet transform acting as detector of the crossing growth singularities. 

• The most relevant quantity provided by the wavelet analysis is the HE of the singularities that we 
determined here through the illustrative example of three microbial species contributing to the fer- 
mentation. 

• In order to have a real technological application of the wavelet techniques in this area, it is important 
to perform calibrations of the mixed-culture growth in terms of their HEs. 

• Elimination of possible numerical artifacts can be achieved by employing wavelets of different vanishing 
moments for analyzing the same set of data. For this goal, we used the Gaussian derivative wavelets, 
which have the same number of vanishing moments as the order of the derivative. It is known from the 
wavelet theory that working with analyzing wavelets with p vanishing moments provide information 
on the regularity of the {p— l)th derivative of the analyzed signal. Since the majority of the numerical 
values of the HEs that we obtained are between one and two, it is clear that the fermentation signal 
under investigation has singularities in the first derivative. At the same time, this explains why we 
did not employed the first Gaussian derivative wavelet since it has less vanishing moments than the 
numerical value of the HEs. 

• It would be very interesting to use other recent techniques to measure the fractality of the total biomass 
signals. We mention the DFA (Detrended Fluctuation Analysis) [51 and the MF-DFA (Multifractal- 
DFA) [6]. These methods are based on determining the scaling properties of the fluctuating moment 
of order q {q = 2 holds for the case of DFA), as a function of the size of the samples and when the data 
contain noises and trends that are unknown from the point of view of their origin and their specific 
forms. In addition, one can consider biological growths whose parameters depend on specific physical 
parameters, such as the temperature, which is common for long-term bacterial activities [7: . Once the 
signal depends on the temperature, the HE depends as well. It would be highly interesting to examine 
the temperature dependence of the HE as implied by the temperature dependence of the total biomass 
signal. 
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